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We present a detailed experimental study of the Uhrig Dynamical Decoupling (UDD) sequence 
in a variety of noise environments. Our qubit system consists of a crystalline array of ^Be"*" ions 
confined in a Penning trap. We use an electron-spin-flip transition as our qubit manifold and drive 
qubit rotations using a 124 GHz microwave system. We study the effect of the UDD sequence in 
mitigating phase errors and compare against the well known CPMG-style multipulse spin echo as a 
function of pulse number, rotation axis, noise spectrum, and noise strength. Our results agree well 
with theoretical predictions for qubit decoherence in the presence of classical phase noise, accounting 
for the effect of finite-duration vr pulses. Finally, we demonstrate that the Uhrig sequence is more 
robust against systematic over/underrotation and detuning errors than is multipulse spin echo, 
despite the precise prescription for pulse-timing in UDD. 



I. INTRODUCTION 



Quantum information systems will be particularly sus- 
ceptible to single bit errors — in stark contrast to their 
digital classical counterparts — due in part to the con- 
tinuous nature of how information is encoded in a qubit. 
While this property makes quantum computers simi- 
lar in some sense to analog computers, the underly- 
ing principles of quantum mechanics and the theory of 
fault-tolerant Quantum Error Correction (QEC) suggest 
that historical considerations about error accumulation 
in analog machines are not applicable to future quantum 
computers [1] . QEC provides a means to suppress qubit 
errors to arbitrarily low levels, as required by the needs 
of a particular algorithm. In order for QEC to provide 
a useful benefit, however, single qubit error probability 
must be sufficiently small that the overhead required to 
implement the error correction procedure does not in- 
crease the likelihood of a logical error more than the ca- 
pacity of the error correction procedure to reduce the 
logical error probability. Suppression of errors below the 
so-called "Fault- Tolerant" threshold ensures a net gain 
in error reduction by application of QEC protocols. Es- 
timated values for the threshold vary widely and depend 
sensitively on the chosen error model and circuit archi- 
tecture, but are generally in the range 10~^ to 10~^ [T]. 
Reaching error rates this low is a challenge for any tech- 
nology, but the requirements of a quantum information 
processor are further complicated by the desire to sup- 
press errors well below this threshold. In the standard 
description, given a qubit error probability p, and an er- 
ror threshold pth, the logical qubit error rate becomes 
Pthip/pthY foi' ^ levels of QEC concatenation (assum- 
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ing a distance 3 code) [I]. Therefore we see that the 
amount by which p is suppressed below ptu determines 
the net benefit in the logical error rate of using a QEC 
code. 

In addition to errors induced by general operational in- 
fidelity, qubits are subject to both bit-flip {X) and phase- 
flip {Z) logical errors induced by coupling of the qubit to 
its environment. These errors arise from two dominant 
forms of decoherence: longitudinal relaxation, and trans- 
verse dephasing, respectively [T]. For example, a qubit 
whose Bloch vector has been rotated from a pole to the 
equator due to longitudinal relaxation has equal proba- 
bilities of being measured |t) or ||) in the qubit basis, 
commensurate with a 50 % probability of a bit-flip er- 
ror. Qubits may couple to many environmental degrees 
of freedom, making decoherence — and especially de- 
phasing — an insidious source of errors for any quantum 
computer. Considerations such as these have motivated 
a search for methods capable of suppressing qubit deco- 
herence, and hence qubit errors [5]. 

Dynamical decoupling is a promising and powerful tool 
in quantum information science for suppressing qubit de- 
coherence [SlIlllSlEllIllHlElllOlin]. This method involves 
the repetitive application of control pulses whose aggre- 
gate effect is to time-reverse the impact of qubit coupling 
to environmental noise. While attractive in its simplicity 
and predicted effectiveness, the use of dynamical decou- 
pling has to date largely been impractical, as most ex- 
perimental quantum information systems have suffered 
from poor operational fidelities in addition to the effects 
of qubit decoherence. However, over the last several years 
trapped ion qubits have demonstrated exceedingly high 
gate fidelities, leading to possible experimental scenarios 
in which qubit errors are dominated by environmental 
coupling rather than poor control. The incorporation 
of dynamical decoupling techniques into a quantum in- 
formation context has been theoretically proposed as a 
means to suppress qubit errors below the fault-tolerance 
threshold, but has been, at this point, only sporadically 
employed for the reasons cited above. 

Understanding of how decoupling sequences suppress 
decoherence is largely based on decades of work in Nu- 
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clear Magnetic Resonance (NMR) [HI US]- It is well 
known from this community that the standard multi- 
pulse spin echo (CPMG, after Carr, Purcell, Meiboom, 
and Gill) technique is capable of suppressing decoherence 
associated with slowly fluctuating noise processes. This 
is, however, an overly restrictive case in a general quan- 
tum informatic context, as environmental coupling may 
contribute noise which fluctuates rapidly compared with 
experimentally relevant timescales such as required qubit 
memory times. Recent theoretical work has introduced 
new methods for understanding the influence of arbitrary 
dynamical decoupling pulse sequences [T31 [THl HSl ttZ] un- 
der realistic experimental conditions, including the use of 
nonidealized control pulses [18], and the presence of noise 
environments with arbitrary power spectra. These tech- 
niques go far beyond the common Magnus expansions |12| 
employed in NMR-style analyses of the influence of en- 
vironmental noise on pulse-sequence performance. How- 
ever, they require detailed experimental veriflcation in 
order to evaluate their utility relative to extant analyti- 
cal techniques. 

Based on these new insights, theoretical [131 [TBI [T7] 
and experimental |18j studies have demonstrated that 
careful manipulation of the pulse-spacing in a dynamical 
decoupling pulse sequence allows improved suppression 
of random phase accumulation over standard, periodic 
pulse sequences [T5J [T3J [T5] in high-frequency-dominated 
noise environments. In particular, Uhrig developed an 
analytically optimized pulse sequence, UDD, which was 
originally designed for qubits subject to an Ohmic noise 
power spectrum in a spin-boson model [TJ. This model 
is of interest to a variety of qubit technologies, includ- 
ing electron spins in quantum dots subject to phonon- 
induced decoherence. Given n 7r-pulses (coherent op- 
erations rotating the Bloch vector halfway around the 
Bloch sphere relative to a defined axis) and a total se- 
quence length r, Uhrig derived an analytical expression 
for the relative positions of the tt pulses such that the se- 
quence strongly suppressed the dephasing effects of noise, 
with analytical and numerical calculations suggesting the 
sequence could outperform the well-established CPMG 
multipulse spin echo |12j by orders of magnitude in the 
high-fidelity (low error) regime. It was later shown that 
the UDD sequence could be employed to suppress deco- 
herence in any noise model |20j . with impressive gains 
resulting in the context of a spin-bath [T^. Together, 
these results suggest that the use of appropriately engi- 
neered dynamical decoupling pulse sequences could pro- 
vide a route to reaching and exceeding the fault-tolerance 
threshold for memory operations without the need to re- 
strict our analysis to slowly-fluctuating noise processes. 

In order to understand the applicability of dynami- 
cal decoupling techniques in future quantum informa- 
tion experiments, we are thus motivated to perform an 
experimental study of the efficacy of dynamical decou- 
pling pulse sequences in the suppression of decoherence- 
induced errors. This study will aim to evaluate theo- 
retical constructs predicting qubit coherence under the 



inffuence of dynamical decoupling pulse sequences, and 
the recent development of novel pulse sequences which 
promise to provide significant benefits under specific en- 
vironmental noise conditions. 

In this article we present detailed experimental stud- 
ies of the performance of dynamical decoupling pulse se- 
quences using an array of trapped-ion qubits as a model 
quantum memory 18, 21,. We drive a ~124 GHz qubit 
transition in a crystallized nonneutral plasma of ^Be"*" 
ions using a quasi-optical microwave system. The ions 
are subjected to environmental noise derived either from 
intrinsic fluctuations in our magnet system, or artificially 
synthesized noise in our control system. We focus on a 
scenario where dephasing dominates longitudinal relax- 
ation; in this framework, appropriate for most contempo- 
rary qubit technologies, dynamical decoupling pulse se- 
quences reduce the accumulation of uncontrolled random 
phases during qubit free-precession periods. Qubit deco- 
herence is measured under application of various pulse 
sequences as a function of sequence length, 7r-pulse num- 
ber n, rotation axis, noise spectrum, and noise strength, 
and is compared against theoretical predictions incorpo- 
rating experimentally relevant parameters. We study the 
performance of the novel UDD sequence and compare 
against the benchmark multipulse spin echo sequence 
(of which CPMG is a special case). Finally, we explore 
the robustness of these pulse sequences in the presence 
of common systematic control errors [22]. Our results 
demonstrate an ability to accurately predict qubit coher- 
ence in arbitrary noise environments under the applica- 
tion of arbitrary pulse sequences, and validate theoretical 
conjectures that careful pulse-sequence engineering can 
indeed provide significant benefits in error suppression. 

We mimic the dephasing processes that may be en- 
countered in solid-state devices by engineering the noise 
environment in our experiment [18j , but gain high-fidelity 
control operations, state initialization, and measurement 
through the use of an atomic system. As such, our exper- 
imental platform behaves as a model quantum system ca- 
pable of emulating the noise environments of other qubit 
technologies in the limit of classical dephasing. This ob- 
servation suggests that our results are generic across all 
qubit realizations and that dynamical decoupling tech- 
niques should yield benefits in all qubit systems. 

This paper is organized as follows. In Section [IT] we de- 
scribe the effect of dynamical decoupling pulse sequences 
on the time-evolution of a qubit subject to random phase 
noise, introducing the concept of the filter function as 
originally derived in references [THl US]. In this section 
we also introduce the construction of the main sequences 
under test, and describe their predicted performance in 
various noise environments. Section IIIII describes our 
experimental system in detail, including discussion of 
qubit structure, coherent control techniques, and hard- 
ware control. Experimental results are presented in Sec- 
tion |IV| comparing the UDD and CPMG pulse sequences 
in 1 /uT and Ohmic noise power spectra, and demonstrat- 
ing strong agreement with theoretical fitting functions. 
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Section |V] studies the accumulation of coherent rotation 
errors due to systematic control infidelities in the appli- 
cation of pulse sequences. Finally, a discussion of experi- 
mental results is presented in Section [Vl] and concluding 
remarks are provided in Section VII 



II. DECOHERENCE UNDER 
PULSE-SEQUENCE APPLICATION 

In order to describe the effects of dynamical decoupling 
pulse sequences on qubit dephasing, we consider classical 
phase randomization due to environmental noise. This is 
described by a Hamiltonian written as 



(1) 



where U, is the unperturbed qubit splitting, /? is a time- 
dependent classical random variable il6j , and az is 
a Pauli operator. Under the influence of this Hamil- 
tonian, a superposition state initially oriented along Y 

evolves in time as \^{t)) = -i=(e~*"*/2g-5 |t) + 

gint/Sgs/o ID) in the lab frame, with |t) and |i) 
the qubit states. In the frame rotating at fi, the term 
I3{t) randomly rotates the Bloch vector in the equato- 
rial plane, corresponding to the accumulation of a ran- 
dom phase between the qubit basis states. The presence 
of any uncontrolled phase increases the probability of a 
logical Z error, and when the root-mean-squared phase 
accumulation is ~ tt the qubit is said to have dephased. 
The characteristic timescale for this process is known as 
T0, and in the absence of relaxation is equivalent to T2, 
generically known as the decoherence time. This form of 
decoherence is known as homogeneous dephasing, as it 
treats the effect of phase randomization on a single qubit 
or in an ensemble where all qubits experience the same 
phase evolution. Alternatively, inhomogeneous broaden- 
ing occurs when qubits in an ensemble are subject to 
different phase evolution, or time-ensemble average mea- 
surements are made, in either case yielding a character- 
istic decoherence time known as r| < T2 [12]. 



A. The Filter Function Formalism 

Hahn showed that in NMR systems injection of a 
TT pulse around X (henceforth denoted -kx) halfway 
through a free-precession period time-reversed the ef- 
fects of inhomogeneous dephasing [23j. Similarly, ran- 
dom phase accumulation can be time-reversed in a ho- 
mogeneous system, so long as temporal fluctuations in 
the noise environment are slow compared to the free- 
precession period. This technique, known as spin-echo 
[I2I [23] . forms the basis of modern studies of dynami- 
cal decoupling |13j . in which multiple qubit rotations are 
used to mitigate random phase accumulation, and hence 
extend qubit coherence time. 



Uhrig [TS] and Cywinski et al. [TB] studied the coher- 
ence of a qubit under the influence of pulse-sequence ap- 
plication in an arbitrary noise environment, and showed 
that for any n-pulse sequence one may write a time- 
domain filter function, 2/„(t), with values ±1, alternat- 
ing between each interpulse free-precession period. The 
time-domain filter function accounts for phase accumu- 
lation under the influence of external noise, with the 
accumulation being effectively time-reversed by applica- 
tion of TTx pulses. Moving to the frequency domain one 
may then write F{lut) — |y„(a;T)p, where yn{u!T) is the 
Fourier transform of the time-domain filter function, and 
T is the total sequence length. 

We modify the time-domain filter function to account 
for finite-length tt pulses [TSJ [M] by incorporating 
a delay with length t-^ and value zero between free- 
precession periods. This approximation assumes that de- 
phasing is negligible during the application of a tt pulse, 
consistent with many experimental observations. Incor- 
porating this delay results in a filter function for an ar- 
bitrary n-pulse sequence 

F{lot) = |y„(c^T)|2 

n 

= |l + (-l)"+ie'"^^ + 2 X;(-l)^e*''^"^cos (c^t^/2)P 

(2) 

where 5jT is the time of the center of the j'^ ttx pulse, 
and T is the sum of the total free-precession time and tt- 
pulse times. The only modification relative to previous 
work associated with accounting for finite r^r comes in 
the addition of the cosine term at the end of the expres- 
sion . The derivation of this equation appears in the 
Appendix of this manuscript. 

The filter function, F{ijjt) provides all information 
about how an arbitrary pulse sequence will suppress 
phase accumulation. Following references [14] [16], for 
arbitrary noise power spectrum Sp{ijj), we write the co- 
herence of a state initially oriented a long Y in a frame 
rotating with frequency O as W(t) = \ {(Jy){'t)\ = e^^^'^\ 
where angled brackets indicate a quantum- mechanical ex- 
pectation value, the overline indicates an ensemble aver- 
age, and 



X{r) 



-^-^ — t (UJTjdcj. 



(3) 



Since the filter function enters the coherence integral as 
a multiplicative factor of S'/3(w), small values of F(ijjt) 
will lead to small values of x(t), and hence coherence 
W{t) - 1. 

This construction describes the infiucnce of an ar- 
bitrary dynamical decoupling pulse sequence on qubit 
coherence in an arbitrary noise environment, and it 
incorporates experimentally realistic conditions such as 
the use of noninstantaneous control pulses. As such 
it represents a significant departure from the Magnus 
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expansion which has historically been employed in 
NMR systems to describe spin coherence in a noisy 
environment. 



B. The CPMG and UDD Sequences 

The multipulse spin echo sequence was originally de- 
veloped in the context of NMR systems where inhomo- 
geneous broadening required an ability to refocus the en- 
semble Bloch vector as it spread out during free preces- 
sion [T2] . However, the same sequence is quite effective 
at suppressing phase randomization due to homogeneous 
effects (e.g. (3{t)) when noise processes are dominated 
by low-frequency components (e.g. Sp{uj) oc I/^jj). An 
n-pulse spin echo sequence has all pulses evenly spaced, 
with the first and last free precession periods half as long 
as the interpulse free precession periods (Fig.[2|:). In the 
case where \'i'o) (the state vector at the leading edge of 
the first 7r-pulse) is oriented along the direction of the ap- 
plied microwave field, the sequence is known as CPMG, 
after its proponents. For reasons that will be described 
in a later section, this particular version of the multipulse 
spin echo sequence is most robust against control error 
accumulation. 

Uhrig showed that manipulation of the relative pulse 
locations, 6j, for fixed n and r leads to modification of 
F(ijjt), providing the ability to tailor the filter function in 
order to strongly suppress noise with certain spectra. In 
particular, he analytically derived an n-pulse sequence in 
which the first n derivatives of yn{^T) vanish for ujt — 0. 
The resulting sequence, UDD, has tt pulse locations de- 
termined analytically as 5j = sin^[7rj'/(2n -I- 2)], for an 
n-pulse sequence (Fig. |2|:). 

The construction presented above allowed Uhrig to tai- 
lor the filter function such that it provided strong sup- 
pression of phase accumulation when noise environments 
possessed significant high-frequency contributions — a 
dramatic advance over CPMG. Uhrig specifically showed 
[15j that in noise spectra including high-frequency com- 
ponents and a sharp high-frequency cutoff, loc, such as 
the Ohmic spectrum {Sp{ijj) oc LoQ{ix) — ujc)) appropri- 
ate for a spin-boson model, the UDD sequence was pre- 
dicted to yield significant gains in performance relative 
to CPMG. By contrast, in the presence of noise with a 
soft high-frequency cutoff, the error-suppression benefits 
arising from the form of the filter function for UDD were 
reduced. These theoretical results were first validated 
experimentally in reference jl8j . 

Numerical simulations showed that by applying the 
UDD sequence, qubit phase error rates for free-precession 
time t, expressed as 1 — W{t), could be orders of mag- 
nitude lower than those achievable using the multipulse 
CPMG sequence when qubits were subjected to high- 
frequency-dominated noise. Accordingly, in this study, 
we will seek to experimentally characterize the perfor- 
mance of the CPMG and UDD sequences in the presence 



of low-frequency and high-frequency dominated noise. 

III. EXPERIMENTAL SYSTEM 
A. Qubit System 

The qubit employed in our experiments is a ground- 
state electron-spin-flip transition of ^Be+ ions in the 
2s manifold, (|t) = |ni/ = 3/2, mj = 1/2) ^ 

Im/ = 3/2, mj = -1/2) = |i)) [18^. This transition is 
highly sensitive to fluctuations in the external magnetic 
field (~28 MHz/mT), but has a negligible spontaneous 
emission rate, making it efficient for studying pure de- 
phasing processes. We form an ordered qubit array of 
~1000 ^Be"*" ions in a Penning ion trap f2T, [25] oper- 
ated at 4.5 T. The motion of the ions along the axis 
of the magnetic field is Doppler cooled to temperatures 
of order 1 mK [27] using ~313-nm UV laser light 
red-detuned from an atomic transition between the 2s 
^5*1/2 and 2p manifolds. When cooled these ions 

form 2D or 3D arrays with well defined crystal structure 
p5l l26l [27l 128] and ion spacing ^10 /im. The ion array 
rotates rigidly at a frequency ^--^20 kHz, due to the ge- 
ometry of the confining electric and magnetic fields. The 
rotational frequency is controlled by an external rotating 
dipole potential |29j [30] . Adjusting the rotation rate of 
the array allows us to controUably manipulate the dimen- 
sionality of the ion crystal (planar vs. 3D). A schematic 
of the experimental system is displayed in Fig. 

B. Initialization, Measurement, and Coherent 
Control 

Qubit state initialization into |t) occurs via optical 
pumping by cooling and repump lasers (Fig. [T^), and 
state readout is achieved by fluorescence detection on 
the same cycling transition used for cooling [31i . In our 
experiments, jf) is a bright state that resonantly scatters 
~313 nm laser light. Fluorescence events are detected 
using a low-dark-count phototube. In any experimen- 
tal cycle, after qubit initialization, the cooling/detection 
laser is pulsed on for a fixed time of 50 ms, fluorescence 
events are detected, a coherent control pulse sequence is 
applied, which under conditions of full coherence results 
in a dark state |J,), and fluorescence is again detected 
for 50 ms. In order to account for optical pumping from 
II) It) induced by the detection laser, the second de- 
tection cycle is broken into five 10 ms bins, and a linear 
fit to the number of counts per bin allows extrapolation 
to a t = count rate (relative to the start of the detec- 
tion cycle). This extrapolated count rate is normalized 
to pre-control-pulse measurements in order to account for 
slow fluctuations in laser intensity. Using these methods 
we readily achieve preparation and measurement of the 
dark-state with fidelity greater than 99.5 %, limited by 
stray-light scatter from the trap electrodes and imperfect 
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FIG. 1: (Color online) Schematic diagram of experimental 
system, (a) Qubit level diagram showing relevant transitions 
and energy levels (see text). (Inset a) Microwave resonance 
transition showing Rabi lineshape for square excitation pulses. 
Vertical axis is normalized counts, (b) Rabi oscillations for 
directly driven qubit transition showing ~ 185 fis. (c) 
Effective Larmor precession measured using a Ramsey in- 
terferometer. Open markers correspond to DDS detuning 
during free-precession time of 100 kHz, solid markers cor- 
respond to 50 kHz detuning, (d) Schematic of microwave 
generation system. CPL=coupler. (e) Schematic of system 
for generating microwave frequency noise, (f) Schematic of 
trap system showing relative positions of hardware elements. 
AOM=Acousto-Optic Modulator. Inset (f) Top-view optical 
micrograph of an ion array showing hexagonal crystal order, 
obtained by imaging ion fluorescence parallel to the magnetic 
field axis, (g) Schematic of experimental control system. 



TTx rotations due to magnetic field inhomogeneities. The 
best measured fidelity of the dark state prepared by a 
single TTx rotation in our system was 99.85 % (optimized 
to maximize microwave coupling to the ions, and mini- 
mize inhomogeneities and stray light scatter). 

Qubit rotations are achieved by directly driving the 
~124 GHz electron spin-flip transition using a quasi- 
optical microwave system. The use of this technique 
eliminates the effects of spontaneous emission which are 
a limiting parameter in Raman-laser-mediated qubit ro- 
tations. Coherent control over this qubit transition 
was first demonstrated in reference [TH]. We gener- 
ate 124 GHz microwaves using a Gunn diode oscillator 
phase-locked to the eighth harmonic of a 15.5 GHz di- 
electric resonator oscillator (DRO), plus a ~77 MHz IF 
reference frequency. The DRO is, in turn, referenced 
to an ultra-stable 100 MHz quartz oscillator (Fig. [l|i). 
The ~77 MHz reference is typically generated using 
a computer-controlled direct digital synthesizer (DDS) 
which may be rapidly switched between four preset out- 
put profiles (each with fixed phase and frequency) via 
TTL inputs. 

Microwaves are launched from a horn antenna and 
propagate through free space before being focused on 
the ion array by a Teflon lens (Fig. Pulsed control 
is afforded via a TTL-controUed p-i-n diode microwave 
switch that provides ~25 dB of attenuation in the off- 
state. The microwave frequency may be detuned while 
the switch is in the off-state in order to minimize the 
effects of microwave leakage. A variable attenuator con- 
trols the amount of microwave power output by the horn 
antenna, and mitigates standing wave reflections from 
the microwave switch that can disrupt the PLL. Ampli- 
tude stability of our microwave system has been mea- 
sured to be >99.9 %. 

A fixed-duration microwave pulse is applied to the 
qubits, and the ~77 MHz reference is swept in order to 
find the resonant transition frequency, /q (Inset, Fig.[ljt). 
This frequency drifts slowly due to flux leakage from our 
superconducting magnet at a relative rate of< 1 x 10~* 
per day. Additionally, /q can occasionally be seen to shift 
suddenly throughout the day by up to ~1 kHz, presum- 
ably due to the influence of external laboratory equip- 
ment in the building. Careful monitoring of /o allows 
us to maintain long-term microwave frequency stability 
relative to the qubit transition of ^1x10^^ over any ex- 
perimental measurement run. 

Varying the length of an on-resonance microwave pulse 
shows Rabi oscillations that decay with a characteristic 
time of 30 to 40 ms (Fig.[T]3). Long-time oscillations are 
best flt with a Gaussian decay envelope arising from mi- 
crowave amplitude instability, but the complex observed 
decay pattern suggests the possibility of several compet- 
ing influences. Tuning the attenuator upstream from the 
horn antenna allows us to vary our ttx time, ~ 70—700 
/iS, while still maintaining reliable ttx pulses. Experi- 
ments in this study were conducted with r^r ~ 185 fis. 
Pulse timing resolution of 50 ns corresponds to a system- 
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atic pulse timing infidelity of < 4 x 10"'''. A tt rotation is 
sensitive to such inaccuracies only to second order. The 
infidelity of vr rotations is ~1 x 10^^, due to magnetic field 
and microwave beam inhomogeneities over the ion crys- 
tal. 

Coherent control of qubit phase is demonstrated us- 
ing a Ramsey interferometer consisting of two {'k/2)x 
pulses spaced by a free-precession period. Microwaves 
are detuned during the interpulse period, resulting in an 
effective controlled Larmor precession, as shown in Fig. 
[TJ; for detunings of 50 and 100 kHz. Ramsey interference 
fringes measured in this way show a characteristic decay 
time of 1 to 2 ms, depending on the size and configura- 
tion of the qubit array. Ion arrays with multiple planes 
generally yield shorter decay times due to axial magnetic 
field inhomogeneities. 



C. Noise Generation 

We synthesize noise in our experimental system in or- 
der to simulate the noise environments of other qubit 
technologies, and to provide a means to validate the 
theoretical constructions employed by Uhrig and others 
[l4l [T5| [TBI lEI • Noise synthesis commences with the gen- 
eration of a desired noise spectrum; all frequency compo- 
nents over a desired range are provided a uniformly dis- 
tributed random phase [0, 27r], and weighted by an enve- 
lope function appropriate for the desired noise spectrum 
(e.g., l/cij). This spectrum is then Fourier transformed, 
yielding a numerical time trace whose ensemble-averaged 
two-time correlation function reproduces the spectrum 
of interest, and sent to an arbitrary waveform genera- 
tor (ARB). The ARB's time-domain voltage output pro- 
vides a frequency modulation signal for a synthesizer, 
whose output is mixed with a DDS (Synthesizer: 275 
MHz, DDS: -197 MHz) to yield the -77 MHz reference 
to the Gunn diode PLL (Fig. ^). 

Our noise synthesis technique involves modulation of 
the microwave driving frequency rather than injection of 
magnetic field noise to the magnet bore. However, modu- 
lation of the microwave frequency used for coherent con- 
trol is equivalent to modulation of the qubit splitting, fi, 
relative to a fixed microwave drive, as either technique 
leads to the same relative phase accumulation between 
the qubit and the control field. As such our noise syn- 
thesis technique directly mimics the dephasing effects of 
external fluctuations in the magnetic field (i.e. f3{t)). 
Further, this technique faithfully recreates the spectrum 
of interest in our system, while direct magnetic field-noise 
injection suffers from, e.g., frequency-dependent shield- 
ing in the open-bore NMR magnet. The noise spectrum 
is characterized from 10^"^ to 10® Hz via direct measure- 
ment of a 20 MHz beat-note produced as described above 
(Synthesizer: 217 MHz, DDS: 197 MHz) using a commer- 
cial phase-noise detection system. 

The noise waveform generated by the ARB and output 
to the FM port of our synthesizer is blocked by a TTL- 



controUed RF switch during ttx pulse application, and 
is asynchronous with the experimental trigger. As such, 
noise is applied only during the interpulse precession pe- 
riod, allowing us to produce high-fidelity 7rx-rotations 
{ttx fidelity is reduced to —95 % if noise is on while driv- 
ing rotations). This technique simulates systems that 
apply very strong, short Tr^-pulses, similar to the origi- 
nal prescriptions for dynamical decoupling. 

In these experiments we set the strength of the injected 
noise spectrum by adjusting the peak-peak amplitude of 
the voltage time-trace output by the arbitrary waveform 
generator, Vn, where the injected noise power scales as 
V^. A representative, smoothed noise spectrum mea- 
sured for V/v =0.7 V is shown in the inset to Fig. |4^. 

D. Experimental Control 

Experiments are performed under software control by 
a PC (Fig.[l^). The PC communicates with hardware via 
various interfaces including USB, parallel port (LPT) and 
general purpose interface bus (GPIB). The central hard- 
ware component of our setup is a programmable logic 
device (PLD) that applies TTL control pulses to other 
hardware, including a DAQ counter for measurement, a 
microwave switch for pulse application, RF switches con- 
trolling repump and cooling lasers, and the profile switch 
of a DDS. In addition to programs sent to the PLD, 
relevant instructions are simultaneously sent to a data 
acquisition board in order to initialize the counter and 
other outputs, and GPIB commands are sent to periph- 
eral hardware. 

For a given experimental sequence, control instructions 
are sent to the PLD via USB. A user command syn- 
chronized to line frequency triggers the PLD to begin 
the experiment. Upon completion of an experiment data 
are read into software from the data acquisition board, a 
reset signal is sent to the PLD, new instructions down- 
loaded, and the next experiment commences. A complete 
experiment including state initialization, measurement, 
and the application of the pulse sequence itself takes ap- 
proximately 200 ms. A single displayed data point typi- 
cally consists of the averaged results of 20-50 experiments 
performed with identical settings. 



E. Pulse Sequencing 

A schematic pulse sequence for an arbitrary, n-pulse 
dynamical decoupling sequence is shown in Figs. |2^ & b. 
The sequence begins with the application of a {tt/2)x- 
pulse which rotates the Bloch vector to the equatorial 
plane. Microwaves are detuned during the interpulse pe- 
riod; application of this detuning results in effective Lar- 
mor precession and also serves to reduce the Rabi rate, 
mitigating unwanted rotations due to microwave leakage. 
Short delays called "Profile Heads (Tails)" between the 
trailing (leading) edge of a pulse and the application (re- 
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moval) of the detuning are employed to precisely control 
the phase evolution of the qubits. For instance, with 100 
kHz detuning, and the first interpulse precession period 
set to an integer multiple of 10 /^s, the Primary Tail may 
be set to 12.5 /iS such that the Bloch vector rotates from 
the Y to the X axis. Varying the length of the Primary 
Tail therefore allows a form of process tomography in the 
equatorial plane to be performed, as will be discussed in 
Section [V] The Final Tail is set such that in the absence 
of dephasing the qubits will be rotated to the dark state 
at the conclusion of the sequence. 

When applying CPMG we select the total free- 
precession time such that each interpulse precession pe- 
riod is an integer multiple of the detuned rotation period, 
10 /is = (100 kHz)~^. Accordingly, in the absence of de- 
phasing the Bloch vector makes an integer number of 
complete rotations about Z during each interpulse pre- 
cession period. Application of UDD generally places no 
constraints on the length of the interpulse precession pe- 
riod, and at the conclusion of each precession period the 
Bloch vector ends at a different location on the equator. 
However, we find that rounding each precession period to 
an integer multiple of 10 ^s such that the Bloch vector 
returns to the same point on the equator before each ttx 
pulse makes no difference in our results within measure- 
ment uncertainty. 

We apply pulse sequences with a few to over 1000 ttx- 
pulses. For long sequences we find that rotation errors 
likely associated with long-time microwave amplitude in- 
stability or magnetic field fluctuations build up, yielding 
minimum measured error rates in dark state probability 
for 1000 TT pulses of ~10 to 15 %. We avoid these errors 
by focusing primarily on sequences consisting of n < 12 
TTx-pulses. In this regime the minimum error rate for any 
sequence is <■ 1 %. 



IV. SEQUENCE PERFORMANCE IN VARIOUS 
NOISE ENVIRONMENTS 



We begin our study by applying the CPMG and UDD 
pulse sequences under the influence of ambient magnetic 
field noise. Ambient field noise is measured directly us- 
ing a field-coil embedded in the bore of our supercon- 
ducting magnet and a low-frequency spectrum analyzer 
(Fig. [3^). We find magnetic field fluctuations varying 
approximately as 1/w^ (^^(w) oc l/ui^), with additional 
sharp spurs, including a prominent feature at ~153 Hz. 
While the exact cause of these spurs is undetermined, 
we believe they are linked to vibrations of the magnet 
system due to both external perturbations and possibly 
liquid nitrogen boil-off in the magnet shroud. We observe 
the ambient field spectrum to vary slowly and discontinu- 
ously in real time, including the overall background level 
as well as the precise positions and heights of the noise 
spurs. The overall scale of the fluctuations in Fig. |3^ 
is calibrated by fitting to the decay of Ramsey fringes 
for a small planar array where effects of inhomogeneous 
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FIG. 2: (Color online) (a) Schematic pulse sequence showing 
relevant timescales. Capitalized letters correspond to specific 
times and relate to images in lower panel. Primary and Fi- 
nal Tails set as indicated in the text. For standard CPMG 
and UDD Primary Tail: 12.5 fis, Final Tail: 17.5 ^s. (b) 
Bloch sphere representation of key experimental procedures 
and events in application of a dynamical decoupling pulse se- 
quence. Curved arrow indicates rotation axis during driving 
or free precession. Driven rotations always about X, free- 
precession always about Z. (c) CPMG (closed markers) and 
UDD (open markers) fractional pulse locations as a function 
of pulse number, illustrating increasing divergence between 
the sequences with larger n. 



broadening can be neglected. 

We apply the CPMG and UDD sequences under ambi- 
ent noise while varying the pulse number, showing that in 
both cases we are able to extend the coherence time with 
n, as predicted by theory (Fig. [3|. Qubit decoherence is 
represented as a colorscale taking values of to 0.5, and 
corresponds to ^(1 — W{t)). Fully dephased qubits have 
an equal probability of being measured in either of the 
two basis states, corresponding to a measured value of 
0.5 for error probability. 

In our experiments we perform tt rotations about the 
X-axis, or produce an effective rotation about Y by com- 
bining a TTx pulse with a net ttz rotation during the free- 
precession period. The latter technique is similar to the 
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XZXZ decoupling pulse sequences designed to suppress 
both longitudinal relaxation and transverse dephasing, 
and may also provide benefits in the context of noise pro- 
cesses that commute with ax- We see that for all times 
and n, these two techniques yield nearly identical results 
despite the need to round the interpulse precession time 
to integer multiples of 10 /is in the case of UDD with Try 
rotations. This similarity suggests that at the 10~^ er- 
ror level (the limit of our combined 7r-rotation and mea- 
surement fidelity), the noise suppressing capabilities of 
UDD are reasonably robust to minor deviations in pulse 
spacing from analytically derived values. Numerical sim- 
ulations confirm that for this noise spectrum, interpulse 
precession times may be rounded by up to 20 ^s before 
the minimum error rate consistently exceeds 10~^. 

While the ambient noise spectrum is convenient as it 
appears naturally in our system, it also provides a useful 
example of a noise power spectrum dominated by low- 
frequency components. In this limit, use of the CPMG 
sequence is predicted to be an efficient means of phase- 
error suppression, and the UDD sequence construction 
should provide no benefits. As expected, UDD appears 
to perform similarly to CPMG in the presence of this 
noise spectrum which has a soft high-cj cutoff. Perfor- 
mance in the high-fidelity (low-error) regime is similar 
between the two sequences, consistent with theoretical 
expressions for qubit decoherence when accounting for 
finite t^, and limits imposed by 7r-rotation and measure- 
ment fidelity. We note, however, that applying UDD with 
odd values of n consistently results in a relatively large 
error rate for times short compared to the coherence time 
for that sequence (~ 10 % vs. ~1 %). This appears for 
both rotation axes and is clearly persistent out to at least 
n = 10, but is not predicted by theoretical expressions 
for qubit decoherence (Figs. [3^ & c). At this point, the 
observed even-odd asymmetry remains unexplained. 

Panels e & f of Fig. |3] show fits to theoretical expres- 
sions for 5 (1 — W(t)) using two parameters: the over- 
all noise strength (a) modifying Sp{Lo) and entering the 
integral x(0 ™ Eq. 3, and the relative strength of the 
153 Hz spur in the measured noise (7). For these fits we 
hold 7 fixed with value 0.15. Theoretical predictions bear 
strong resemblance to the experimental data incorporat- 
ing the measured ambient Sfj{Lo). As such, fitting factors 
of order unity are expected. In the absence of significant 
TT-pulse imperfections, the same theory applies for -kx 
and TTy rotations. Fits are performed for each value of n 
and show minor variations in extracted a between traces. 
The mean value of the fitting factor is a =2.3±0.3, with 
individual fit uncertainties in the range of 0.1-0.3. Vari- 
ations in extracted values of a are consistent with real- 
time observation of changes in the noise spectrum. We 
believe that deviations between our experimental data 
and fitting functions are dominated by slow as well as 
discontinuous changes in the ambient noise environment. 

At the other extreme of possible noise power spectra, 
we consider noise dominated by high-frequency contribu- 
tions for which the simple CPMG sequence is not known 
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FIG. 3: (Color online) (a-d) Colorscale images of measured 
qubit decoherence under ambient noise as a function of to- 
tal free precession time and pulsenumber, n, for various se- 
quences. Colorscale illustrates phase error probability as nor- 
malized counts (0.5 = 50 % probability of ending sequence as 
IT) or ID), with light colors corresponding to low qubit errors 
and black corresponding to full dephasing. Axis of rotation 
is designated in the figure. Try rotations driven by a combi- 
nation of a TTx pulse and a net ttz rotation during each free 
precession period, (e-f) Theoretical fits to qubit decoherence 
using the filter function as defined in Eq. (2), and the noise 
spectrum plotted in panel (g). (g) Measured noise spectrum 
showing \/uj^ field fluctuations (5/3(01) oc l/<^*). Note strong 
spur at ci; ^ 2Tr X 153 Hz. Dashed line is a guide to the eye 
representing S0{ui) oc l/o;* scaling. 



to be effective. We mimic the spectrum produced by a 
spin-boson model in the classical limit for which UDD 
was originally developed by synthesizing classical noise 
with an Ohmic spectrum and a sharp-high frequency cut- 
off [32]. Numerical simulations suggest that with this 
noise spectrum, in the high-fidelity regime UDD is capa- 
ble of significantly outperforming CPMG, yielding rela- 
tive gains of several orders of magnitude [14] I16j. 

In our experimental system the minimum observable 
error rates in the high-fidelity regime are limited by 
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measurement infidelity ^^0.5 %, making differences in 
sequence performance indistinguishable in the regime 
Uhrig originally studied. In order to overcome such lim- 
itations we emphasize the differences between these se- 
quences by synthesizing noise with high power relative 
to the ambient spectrum, thus raising the minimum er- 
ror rates to above the ~1 % level. This approach is 
well-suited to a study of UDD as follows. In his orig- 
inal study, Uhrig demonstrated that for fixed resources 
and an Ohmic noise spectrum, significant error suppres- 
sion could be gained in the high-fidelity regime at the 
expense of a slight decrease in coherence time relative to 
other sequences [11]. This is manifested as a "crossover" 
point about which the relative performances of UDD and 
a benchmark sequence (e.g. CPMG) are expected to be- 
come inverted. If this crossover point occurs at low er- 
ror rates in the high-fidelity regime the benefits of UDD 
will, in our system, be swamped by measurement infi- 
delity. Uhrig also showed, however, that increasing noise 
strength can push the crossover point to higher error 
rates, thus allowing us to enter an experimental regime 
where differences in sequence performance are easily dis- 
tinguished. Thus, in our experimental study we utilize an 
artificially engineered Ohmic noise spectrum with noise 
power many orders of magnitude larger than the ambient 
spectrum. While our noise does not contain quantum me- 
chanical correlations, the generality of the optimization 
procedure which led to the discovery of UDD is tested by 
employing such a noise spectrum. 

A smoothed power spectrum with Ohmic scaling, mea- 
sured for Vn — 0.7 V is displayed in the inset to Fig. |4^. 
The measured spectral peak at the cutoff of the ohmic 
spectrum with Vn = 0.7 V is approximately five orders of 
magnitude larger than the ambient spectrum at the same 
frequency. We apply the CPMG and UDD sequences in 
the presence of this noise power spectrum for Vn =2 V 
and display the results for various values of n in Fig. |4^. 
The UDD sequence clearly and dramatically outperforms 
CPMG for all measured values of n. Error bars are shown 
for data with n = 6, and are derived from the statistical 
variations in the individual experiments that are aver- 
aged to produce a single data point [33j . demonstrating 
that the traces differ with statistical significance. As ex- 
pected, the coherence time may be extended with the 
addition of ttx pulses. Data show good agreement with 
theoretical predictions for qubit decoherence using the 
measured noise spectrum and a single fit parameter, a, 
the overall noise scaling. 

We show the performance of UDD and CPMG for var- 
ious noise strengths, parameterized by Vn, in Fig. ^p. 
As predicted in the discussion above, increasing the 
noise power results in a reversal of the relative perfor- 
mance of the two sequences, with UDD's performance 
relative to CPMG improving with Vn- Again, we find 
good agreement between data and theory. In particular, 
we expect an extracted fit parameter close to unity for 
Vn = 0.7, the value of ARB output at which the noise 
power spectrum is measured. We find ajjoD = 0.88 and 
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FIG. 4: (Color online) UDD (open markers) and CPMG 
(solid markers) performance under the application of an 
Ohmic noise spectrum with a sharp high-w cutoff, (a) UDD 
and CPMG for fixed noise power and various values of n. 
Representative error bars presented for n = 6. Fits show good 
agreement with data (see text). Vertical axis corresponds to 
an experimental measure of |(1 — W{ty) as in previous fig- 
ures. Data offset vertically for clarity. Inset (a) smoothed, 
measured ohmic noise spectrum with high-oj cutoff at 27r x 
500 Hz. Actual spectrum used for fitting extends from 27r x 10 
mHz to 27r X 1 MHz, truncated for visibility, (b) Evolution of 
pulse sequences with increasing noise strength, parameterized 
by Vjv , showing that UDD performance relative to CPMG im- 
proves as noise intensity increases. Data offset vertically for 
clarity, (c) Quadratic scaling of extracted fit parameters for 
various noise strengths as a function of Vat on a log-log plot. 
All fit parameters normalized to extracted fit parameter for 
Ujv = 0.1 V. Solid line represents quadratic scaling of the nor- 
malized value of a with Vjv. Data for both CPMG and UDD 
presented. 

ctCPMG = 0.87. The small, consistent deviation of these 
values from one is likely due to the difference in carrier 
frequency employed in phase noise detection versus noise 
injection. 

The validity of our fitting procedure is further vali- 
dated by plotting avjv/ao.i (Fig- the extracted fit 
parameter for each value of Vn normalized to that ex- 
tracted from Vn = 0.1 V, as a function of Vn- Similar- 
ity of the extracted fit parameters for CPMG and UDD 
is expected as data are taken for both sequences under 
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nominally identical conditions. These values lie along 
a straight line with slope two on a log-log plot, con- 
sistent with noise power that scales quadratically with 
Vat (noise power scales as the square of the field fluc- 
tuations). As an example, for Vn — 0.5 V, we expect 
ao 5/ao.i = (0.5 / 0.1)2 ^ 25, and obtain for UDD 30±5. 
The strong agreement between data and theory is indica- 
tive of the general capabilities of the presented theoretical 
framework. 



SEQUENCE ROBUSTNESS TO 
SYSTEMATIC ERRORS 



Systematic errors associated with quantum control 
mechanisms promise to be a significant challenge for re- 
alizing a quantum information system [1, . It is therefore 
imperative to study the robustness of any dynamical de- 
coupling pulse sequence against errors associated with 
the application of imprecise Tr^-pulses [HI HI]. Further, 
as quantum information systems will require that dynam- 
ical decoupling techniques be applied to arbitrary initial 
states, it is vital that the sequence performance be known 
for a variety of input states, information gained through 
a procedure known as process tomography. In this sec- 
tion we present measurements of the effect of rotation 
errors on UDD and CPMG sequence performance, char- 
acterized for a variety of input states on the equatorial 
plane of the Bloch sphere (following the initial {tt/2)x 
pulse) . 

In NMR spectroscopy, the CPMG sequence is known 
to be particularly robust to over- and under-rotation er- 
rors associated with inaccurate ttx pulses [T^]. The se- 
quence, however, is highly sensitive to the initial state 
on the equatorial plane, providing robust performance 
only when the Bloch vector is aligned with the direction 
of the applied microwave field. In this situation, with 
1^0) = |Ar), the application of a ttx pulse yields only 
a small excursion of the Bloch vector about X. By con- 
trast, initializing |^'o) = \ Y), the application of ttx pulses 
produces large excursions of the Bloch vector along the 
prime meridian of the Bloch sphere, and provides little 
resistance to rotation errors. Similarly, due to the nature 
of the pulse spacing in UDD, when the interpulse pre- 
cession periods are unconstrained (and unrounded) we 
find that for each Tr^-pulse the Bloch vector assumes a 
different location on the equator and rotates in a plane 
parallel to that defined by the prime meridian. As such 
the UDD sequence does not immediately appear as if it 
will provide the same robustness to pulse errors as does 
CPMG. 

We perform process tomography for n = 6 ttx pulse 
UDD and multipulse spin echo under the influence of 
ambient noise. The initial state is varied by setting the 
"primary tail" (see Fig. |2^) in our pulse sequence such 
that the Bloch vector rotates by a fixed amount prior to 
the application of the first ttx pulse. A rotation about 
Z by angle 7r/2 corresponds to the Bloch vector initially 
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FIG. 5: (Color online) Robustness of n = 6 spin echo and 
UDD pulse sequences to systematic rotation errors for a va- 
riety of input states, measured for free-precession time 8 ms 
and ambient noise, (a-b) Colorscale representation of qubit 
rotation to ||) as a function of microwave reference frequency 
(centered on /o) and initial state. Dark colors indicate rota- 
tion errors, (c) Difference between UDD and echo results for 
above parameters. White indicates regions where UDD out- 
performs spin echo, (d-e) Colorscale representation of normal- 
ized counts as a function of actual "vr-pulse" rotation angle 
(tt =185 fis) where is manipulated to mimic systematic 
rotation errors, (f) Difference between UDD and echo results 
for "TT-pulse" rotations, white shows regions where UDD out- 
performs echo. Randomly located black dots correspond to 
laser dropouts during data collection, (g) Bloch sphere repre- 
sentations of several relevant starting states used in colorscale 
images above. 



pointing along X, parallel to the direction of the mi- 
crowave field vector. In these experiments we induce ro- 
tation errors by two distinct means: changing the length 
of the TTx pulse relative to the actual r^, and detuning the 
microwave frequency from resonance. These two error 
mechanisms inject systematic over- and under-rotations 
as well as phase rotations. 

Our measurements clearly demonstrate that the UDD 
sequence provides significantly more robustness to errors 
over a wider variety of states than does multipulse spin- 
echo (Fig. [5]). We observe that the CPMG sequence (ini- 
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tial state rotated {1^/2) z) provides the best performance 
under both forms of systematic error, yielding error rates 
less than about 10 % for deviations up to 50 % of r^r or ap- 
proximately half a resonance linewidth. The multipulse 
spin echo sequence, however, is far less robust to error 
for other initial states on the equator. By contrast, the 
UDD sequence shows strong robustness across a variety 
of initial states. For example, for a ^--^14 % deviation in 

and an initial rotation tt^, the UDD sequence shows a 
final state rotation error of '--^lO % while the multipulse 
spin echo has been overrotated by ~70 %. Similarly, a 
systematic detuning of 500 Hz (1/4 linewidth) yields an 
accumulated error of ^8 % for UDD while multipulse 
echo experiences an error of ~50 %. Panels c & f of 
Fig. [5] show the difference between results for the two 
sequences, with white indicating conditions under which 
UDD outperforms multipulse spin echo. 

These results are striking, as they indicate that de- 
spite the precise prescription for the conditions under 
which UDD is shown to work well, it appears to be quite 
robust to systematic errors, and even more robust than 
standard spin echo for a variety of input states. Errors 
may be suppressed to the few-percent level with UDD, 
even for sizeable systematic rotation errors. While such 
errors are not due to a loss of phase coherence, they rep- 
resent large sources of logical phase- flip and bit-flip errors 
in the case where dynamical decoupling pulse sequences 
are incorporated into more complex quantum algorithms. 



VI. DISCUSSION 

A key requirement for the adoption of any quantum 
control technique is the ability to accurately predict a 
qubit's state after the application of that technique. In 
the context of dynamical decoupling, that requirement 
is translated to an ability to accurately predict qubit 
coherence under pulse sequence application. Historically, 
this capability was derived from experimentally vali- 
dated procedures borrowed from NMR, but significant 
limitations prevented general applicability. 

The sensitivity of quantum information processing 
tasks to dephasing-induced errors requires more ac- 
curate techniques for predicting qubit coherence, and 
more general error-suppression strategies than those 
introduced in NMR. In this work we have described a 
general formalism to predict qubit coherence under the 
application of dynamical decoupling pulse sequences 
with noninstantaneous control pulses. This formal- 
ism applies broadly to arbitrarily spaced pulses, and 
accounts for the influence of arbitrary noise spectra. 
Moreover, we have experimentally implemented UDD, a 
novel dynamical decoupling pulse sequence constructed 
using the above formalisms in order to provide enhanced 
error-suppression performance in the presence of fast 
noise. 

Our experimental results validate the conditions 
under which the UDD sequence is predicted to outper- 



form CPMG, and in so doing, validate the underlying 
theoretical constructs predicting qubit coherence. 
We have studied the relative performance of these 
sequences under the approximate limiting cases of 
low-frequency-dominated {Sp{u}) oc 1/^^) and high- 
frequency-dominated {Si3{lo) oc ljQ{ijJ — loc)) noise power 
spectra, in both cases obtaining good agreement between 
measurements and theoretical calculations. We expect 
sequence performance to be intermediate for noise power 
spectra falling between these approximate limiting cases, 
as is likely for experimental systems which suffer from a 
variety of noise processes. 

Our measurements demonstrate that in the presence 
of noise power spectra with significant high-frequency 
contributions, the UDD sequence is capable of outper- 
forming the benchmark multipulse spin echo (CPMG). 
These measurements employ artificially synthesized noise 
with very high power in order to overcome measurement 
limitations, but there is no reason to believe that the 
agreement between theory and experiment would break 
down in the low-error regime in the presence of weaker 
noise, or in a system with improved measurement 
fidelity. Accordingly, we believe that the UDD pulse 
sequence — which requires no more control resources 
than CPMG — may provide benefits of many orders of 
magnitude in the presence of appropriate noise spectra. 

We remind the reader that measurements are limited 
in the high-fidelity regime by readout infidelity. We 
emphasize that the specific observed error rates or coher- 
ence times are not the main results of our study. Instead, 
the precise ability of theory to replicate coherence line- 
shapes for various pulse sequences, pulse numbers, noise 
spectra, and noise strengths demonstrates the general 
applicability of these techniques. 

Even in the case of low-frequency-dominated noise, 
our measurements and previously presented numerical 
simulations [m [151 HSl Hi] suggest that the performance 
of the UDD sequence is similar to that of CPMG (rather 
than substantially worse), except in the measured 
coherence time. For long free-precession times, in low- 
frequency-dominated noise environments CPMG indeed 
outperforms UDD. In a computational setting, however, 
the accumulation of error near the 1/e decay level is 
already many orders of magnitude larger than that 
tolerable in a fault-tolerant context. Accordingly, from 
a practical perspective, UDD may remain an attractive 
choice for an experimentalist with an ill-defined or 
discontinuously variable noise environment, even if the 
noise is typically dominated by low-frequency compo- 
nents. This suggestion is supported by the measured 
robustness of the UDD sequence to coherent rotation 
errors and timing inaccuracies across a variety of relevant 
input states. Further, this robustness may make UDD 
an attractive sequence for NMR spectroscopists. 
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VII. CONCLUSION 

In this article we have provided a detailed experimen- 
tal study of the application of dynamical decoupling 
pulse sequences in a model atomic system. In particular, 
we have compared the performance of an analytically 
optimized sequence, UDD, to standard multipulse 
spin echo in a variety of noise environments. Our 
measurements have verified theoretical predictions 
that under certain experimental conditions, UDD can 
dramatically outperform the familiar CPMG sequence. 
Further, we have been able to accurately predict qubit 
decoherence using an extension of theoretical constructs 
that incorporates nonzero t.^. Fits using our modified 
theoretical expressions show strong agreement with 
experimental data as a function of free-precession time, 
TT pulse number, noise spectrum, and noise strength. 

In addition to verifying a large body of literature 
on the effects of dynamical decoupling generally, and 
the performance of the UDD sequence specifically, 
our experiments solidify the utility of dynamical de- 
couphng as a viable means of error suppression for 
quantum information processing. The high degree of 
agreement observed between data and theory, coupled 
with previous theoretical calculations, suggest that the 
use of optimized dynamical decoupling pulse sequences 
has the ability to suppress qubit phase-flip errors to 



near or below the fault-tolerance threshold without an 
exponential increase in resources. 

Future experiments will explore the benefits of con- 
catenation of optimized pulse sequences in suppressing 
phase errors, and the use of composite pulses to improve 
control fidelity. In addition, we will study more gen- 
eral noise models including the effects of longitudinal 
relaxation and quantum mechanical correlations by the 
coherent excitation of phonon modes in our ion crystal. 



APPENDIX: DERIVATION OF FILTER 
FUNCTION ACCOUNTING FOR 
FINITE-LENGTH tt PULSES 



In this appendix we present the derivation of a modi- 
fied filter function that accounts for the infiuence of finite 
T^, as used in Eq. (1). 

We introduce the following convention: Assume the 
total dynamical decoupling sequence length, delays plus 
the sum of all 7r-pulse durations, is r. If the center of 



the j'th TT-pulse occurs at time t = tj then let Sj = tj/r, 
Sn+i = 1 and let each 7r-pulse have duration 
The filter function is given by the modulus squared of 
y-ni^T) which is defined as 
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0, for t > 5n+lT 



(between 7r-pulses) 
(during 7r-pulses) 



(A.2) 



then 



y„{ajT) = |- 



Juj{Si-S.^/2)t _ 



^iLu{S2-S^/2)T _ ^i0j{Sl+S^/2)T 



= |l-e''^^i'" 



+ e 



+ ... +(-1) 

iU}62T 



■]} 



giaJ(5-7r/2T _|_ ^—icod-n- /2t | 



+ ... +(-l)"+^ne^ 



1 



1 + (_i)(»+i)e*"'5"+i^ + 2j2{-iye"^^^^ cos (wr^/2) 



(A.3) 
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